Asymptotic analysis of the Gunn effect with realistic boundary conditions. 
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A general asymptotic analysis of the Gunn effect in n-GaAs under general boundary conditions 
for metal-semiconductor contacts is presented. Depending on the parameter values in the boundary 
condition of the injecting contact, different types of waves mediate the Gunn effect. The periodic 
current oscillation typical of the Gunn effect may be caused by moving charge-monopole accumu- 
lation or depletion layers, or by low or high-field charge-dipole solitary waves. A new instability 
caused by multiple shedding of (low field) dipole waves is found. In all cases the shape of the cur- 
rent oscillation is described in detail: we show the direct relationship between its major features 
(maxima, minima, plateau's, . . . ) and several critical currents (which depend on the values of the 
contact parameters). Our results open the possibility of measuring contact parameters from the 
analysis of the shape of the current oscillation. 
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I. INTRODUCTION 



The Gunn effect appears in many semiconductor sam- 
ples presenting negative differential resistance and sub- 
ject to voltage bias conditions [0-0. It consists of a 
periodic shedding of pulses of the electric field at the 
injecting contact, which then progress and are annihi- 
lated at the receiving contact. As a result there appears 
a periodic oscillation of the current through a passive 
external circuit attached to the semiconductor. Under 
different conditions, the current self-sustained oscillation 
may be caused by the motion of charge accumulation lay- 
ers (charge monopoles) [p|, not by the usual electric field 
pulses which are charge dipoles. Most of the experiments 
on the Gunn effect in different materials take place in 
samples with attached planar contacts, so that the wave 
motion may be safely assumed to be one-dimensional. 
Despite the vast literature on the Gunn effect, it is sur- 
prising that many basic questions remain poorly under- 
stood. For example, given a description of the charge 
transport in the bulk semiconductor (say at the level of 
drift-diffusion and rate equations), which are the proper 
boundary conditions for given contacts and how they af- 
fect the self-sustained current oscillations. The first ques- 
tion has been addressed in a companion paper, Ref. ||, 
while the second will be answered here. 



Until recently, when confronted with the Gunn effect, 
theorists would resort to computer simulations of more 
or less complicated models (which were supposed to re- 
flect the physics of a given semiconductor), and would 
then explain qualitatively their numerics. Resort to spe- 
cial solutions valid for infinite semiconductors at constant 
current bias conditions || , or to extrapolations of Kroe- 
mer's NL criterion were often used to interpret the 
simulation results. This left the processes of generation 
and annihilation of domains at the contacts (and in fact 
it also left out the dynamics of wavefronts and pulses) 
outside theoretical considerations. Concerning asymp- 
totic descriptions of the Gunn effect which delve deeper 
than just numerical simulations of drift-diffusion models, 
some progress has been made recently [p~0|— p~3| - These 
works propose asymptotic descriptions of the Gunn effect 
exploiting the fact that this effect is seen most clearly 
in semiconductors having a large value of the product 
of sample length times doping (basically a dimensionless 
length). The role of the "N-L product" in the analysis 
of the Gunn instability was already discovered by Kroe- 
mer [p| and exploited to study the linear stability (small 
signal analysis) of stationary solutions by many authors 
ra. It was recognized only much later that in the limit 
of large dimensionless length (N-L product) it is possi- 
ble to describe asymptotically both the onset jl6| and 
the fully developed Gunn instability [ftOl. In this asymp- 
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totic limit the processes of repeatedly generating a new 
wave (a charge monopole or dipole domain) at the inject- 
ing contact, the motion of the wave and its annihilation 
at the receiving contact may be well separated. Then 
they can be analyzed and combined to fully describe the 
Gunn effect. In particular, the effect of contacts on these 
processes and in determining the shape of the current 
oscillation can be clearly stated. In this paper we use 
our asymptotic theory to study Kroemer's model for n- 
GaAs under b.c. corresponding to ideal MS contacts. We 
find that these b.c. give rise to a multivalued control 
current-field characteristic at the injecting contact. The 
asymptotic analysis shows that the Gunn effect can be 
mediated by both charge monopole or dipole domains 
according to the values of the contact parameters. Shed- 
ding of new charge dipole waves from the injecting con- 
tact is adiabatic in clear distinction with what happens if 
the control characteristic of the contact is single-valued 
[Ejyni. In the later case (analyzed for a p-Ge model 
in Ref. ]I^)) the charge dipole pulses are created very 
rapidly at the injecting contact and they advance and 
grow simultaneously (see also |l(J), whereas for multi- 
valued control characteristic the boundary layer at the 
injecting contact grows adiabatically to a much greater 
size before a new pulse can be shed. These facts may de- 
termine appreciably the shape of the current oscillations. 
Our analysis could be extended to more complex mod- 
els also displaying the Gunn effect, []l2l-[l4||. Depending 
on the values of the parameters characterizing the in- 
jecting contact and of the dc voltage bias, we find Gunn 
oscillations mediated by charge accumulation and deple- 
tion monopole wavefronts, and high and low field charge 
dipole domains. We also find narrow regions in the pa- 
rameter space where multiple shedding of dipole domains 
at the injecting contact occurs. We have thus found a 
characterization of all possible dynamic behaviors which 
an ideal metal-semiconductor contact would give rise to 
in the Kroemer model for the Gunn effect. This opens the 
possibility of extracting information about the contacts 
from the analysis of the Gunn oscillations themselves, a 
subject of considerable interest for applied researchers. 

The rest of the paper is as follows. In Section II we 
present Kroemer's model and the boundary conditions 
for metal-semiconductor contacts discussed in the com- 
panion paper [gj. In Section III we present our asymp- 
totic analysis of Kroemer's model and find that different 
types of Gunn effect are possible according to the values 
of the bias and of certain dimensionless parameters ap- 
pearing in the b.c, io and ao- charge monopoles (moving 
charge depletion and accumulation layers), high and low 
field solitary waves (charge dipoles), multiple (low field) 
charge dipoles, are predicted and confirmed by numerical 
simulations. Section V contains our conclusions whereas 
the Appendices are devoted to different technical matters 
related to the main text. 



II. KROEMER'S MODEL AND B.C. FOR 
METAL-SEMICONDUCTOR CONTACTS 



The unipolar drift-diffusion model for the Gunn effect 
proposed by Kroemer d]l||] , is generally accepted to pro- 
vide a rather complete description of the main features 
of this effect. In the dimensionless units described in Ref. 
Ipfl, Kroemer's model is 



dE , . (dE 
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E(x, t) dx — (f>. 



(1) 
(2) 



Eq. (0) is Ampere's law which says that the sum of dis- 
placement current and drift-diffusion current is equal to 
the total current density J(t). It can be obtained by dif- 
ferentiating the Poisson equation, dE/dx = n — I, with 
respect to time, substituting the charge continuity equa- 
tion dn/dt + dj(x, t)/dx = [the electron current density 
is of the drift-diffusion type: j(x,t) = nv(E) — 5dn/dx], 
and then integrating the result with respect to x. The 
electron velocity is assumed to be N-shaped. For specific 
numerical calculations we shall use Kroemer's curve Il31 



v(E) = E 



l + BE i 
l + E A ' 



(3) 



(it has a maximum vm > at Em > followed by a 
minimum < v m < vm at E m > Em) and the electron 
difusivity S to be constant. The dc bias <f> is the average 
electric field on the semiconductor sample. Eqs. (|l])-(H) 
need to be solved with an appropriate initial field pro- 
file E(x, 0) and subject to the corresponding b.c. For 
ideal metal-semiconductor the following mixed boundary 
conditions have been derived in Hi 



dE , s / , % dE , s 
— (0,f) = a li -J(t) + -^r{0,t) 



(4) 



^(L,t)=aJi L + J(t)-^(L,t)) C) 



where ii and cti, (i = 0,£), are dimensionless parame- 
ters which are a combination of the semiconductor effec- 
tive density of states, contact barrier height, Richardson's 
constant, doping and temperature (see ||). In what fol- 
lows, it will be assumed to be positive because the phys- 
ically interesting phenomena (including the usual Gunn 
effect mediated by high field domains) are observed for 
these values of ii (see phase diagram in || ) . 

For typical n-GaAs data, S < I and i » 1 §. In 
this limit, we shall find approximate solutions to the 
initial boundary value problem (0)-((|) for E(x,t) and 
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J(t). Strictly speaking, the simple asymptotic descrip- 
tion that follows holds in the limit L — > oo, even when 
S = 0(1) |l3[| . Assuming i5 < 1 just simplifies the de- 
scription of the traveling waves of electric field in the 
semiconductor through the use of characteristic equa- 
tions and shock waves [|lO|Jl7| , [l^1 . For example, it is 
shown in Appendix ^ (by using the method of character- 
istics) that the boundary condition (^) implies that the 
electric field at the injecting contact, £7(0, t) = Eo(t), 
obeys the following equation: 



dt 



3c(E) = 



= J-j c (E ), 

(1 + a i )v(E) 
1 + a v(E) 



(0) 
(7) 



These expressions constitute a Dirichlet boundary condi- 
tion for the electric field which contains the same infor- 
mation as the mixed condition (|j). The contact curve, 
ic (-£■)> presents two extrema, a minimum, j cm — j c (E m ), 
and a maximum, jcM — 3c(Em), at the same field val- 
ues as the electron velocity curve v(E). j c {E) tends 
to j$ at — a^ 1 + iq for high electric fields. As we 
shall see below, during most of an oscillation period, 
jc{Ea) ~ J, and this expression yields a multivalued 
contact-characteristic curve relating the field at the in- 
jecting contact to the actual value of the current density. 

To take advantage of the large-length limit, we will use 
the following rescaled time and length, 



1 

e = — . 
L 



t 
L 



V= L 



Then Eqs. become 



i dE rvs dE \ ,2^ E 



(8) 



(9) 



E(y,s)dy = 0. (10) 



Notice that the boundary condition ^ becomes 

dE 



J-j c {E)=e- 



ds 



(11) 



in the limit 6 -C 1. 



III. ASYMPTOTICS OF THE GUNN EFFECT 



In a previous paper Q we have analyzed the stationary 
solutions of Kroemer's model with metal-semiconductor 
contacts and discussed their stability. No stable station- 
ary solution is expected for certain ranges of bias and 



the cathode contact parameters io and ao- In these cir- 
cumstances, the Gunn effect mediated by either mov- 
ing charge monopoles or dipoles might appear. A rich 
phenomenology of propagating waves and current oscil- 
lations has been numerically observed for these param- 
eter ranges. Among them, we have observed both high 
(Figs. [j], ||) and low (Fig. ||) field solitary wavesfmoving 
charge dipoles) , multiple low field dipoles (Fig. , mov- 
ing charge accumulation (Fig. ||) and charge depletion 
monopoles (Fig. ||). In Ref. 

citecompanion, we have identified the critical currents 
determining which type of waves mediate the current os- 
cillation (in the limit S -C 1). They are related to the 
boundary conditions in the following way: 

• H jcM > vm 1 the Gunn effect is mediated by mov- 
ing charge accumulation monopoles [j C M and vm 
are the local maxima of the contact current and 
electron velocity curves, both reached at E = Em]- 

• If jcm < Vm i and j'q * > v m , the Gunn effect is me- 
diated by moving charge depletion monopoles [j cm 
and v m are the local minima of the contact cur- 
rent and electron velocity curves, both reached at 



E 



TP . -sat 
-C'mj JO 



io is the value at which 



j c (E) saturates at high electric fields]. 

• If v m < j cm < j cM < v M , moving charge dipoles 
mediate the Gunn effect (see Ref. || for more de- 
tails). 

As shown in the figures, there are several stages in 
each period of the current oscillation corresponding to the 
processes of generation, propagation and annihilation of 
electric field domains. Each stage has its own time and 
space scales, and hence the oscillation can be suitably 
described by means of a matched asymptotic analysis. 
For instance, the annihilation of wavefronts takes place 
on a fast time scale compared to that governing wave- 
front propagation. Thus on the time scale of wavefront 
propagation, the annihilation of wavefronts is a quasi- 
instantaneous process during which the time derivative 
of the current density changes appreciably while the cur- 
rent itself, J(s), does not change. On the other hand, the 
generation of fronts takes place adiabatically on a much 
slower time scale comparable to wavefront propagation. 
In this stage both J(s) and its derivative change appre- 
ciably. This is quite different from the fast generation of 
fronts and pulses observed for other types of boundary 
conditions |12|. 



As long as these different processes take place on differ- 
ent time and space scales, there are different stages of the 
oscillation which can be analyzed separately. This hap- 
pens for certain bias ranges. We shall then construct the 
approximate electric field and current density solutions 
by means of matched asymptotic expansions. A detailed 
description of a period of the current oscillation will be 
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then obtained. For other bias values, several processes 
occurs almost simultaneously (e.g., the annihilation of a 
wavefront may occur during the process of detachment 
of another wavefront from the injecting boundary layer 
for low bias values). This complicates the asymptotic 
description without adding much to our physical under- 
standing, so that we will omit the details. 

Note that in the limit e = 1/L -C 1, the solutions 
of Eqs. (P|)-(|To|) are piecewise constant: on most of the 
y-interval, E is equal to one or another of the zeros of 
v{E) — J [notice that this equation may have three ze- 
roes, we denote by E\ < E 2 < E 3 ], separated by transi- 
tion layers that connect them. At y — and y = 1 there 
are boundary layers (quasi-stationary most of the time), 
which we call injecting and receiving layers, respectively. 
It can be seen that the propagation of fronts turns out to 
be a quasi-stationary process while the generation is not. 
Thus, two different asymptotic approaches will be used: 
one for the description of the quasistationary propagation 
of fronts and the other for their (non-quasistationary) 
generation. 



A. Quasistationary propagation of wavefronts 



In this section we will present the asymptotics of the 
quasistationary propagation of wavefronts, and of the 
corresponding time evolution of the density current. 

Wavefronts are moving transition layers connecting re- 
gions of the sample where the electric field is spatially 
uniform. In order to describe their quasistationary prop- 
agation we will proceed as follows. We assume that 
E{y, s) is either E\(J) or E 3 (J) outside boundary layers 
and wavefronts. Let the wavefront located at y = Y(s) 
move with velocity c = dY/ds. For each value of J, the 
wavefront advances with speed c = c+(J), if it connects 
E = E X {J) [y < Y(s)\ to E 3 (J) [y > Y(s)], or with speed 
c = c_(J), if it connects E = E 3 (J) [y < Y(s)] to E X {J) 
[y > Y(s)] . To find the inner structure of a wavefront and 
its speed, we introduce a coordinate £ = e _1 [y — li(s)] 
moving with the wavefront. Then we need to solve the 
following problem for the equation 



dE -p. 



dF 



[v(E) ~c]F + v(E) - J 



(12) 



obtained from Eq. (^): Find the unique value c — c + (J) 
[resp. c_(J)/ such that there is a solution of |7j| ) with 
E(-oo) = E^J) and E{po) = E 3 (J) [resp, E(-oo) = 
E 3 (J) and E(oo) = E\{J)]. The solution of this prob- 
lem will provide both the speed and inner structure of 
the wave front. In terms of the phase plane (|l2|), the 
previous problem is equivalent to find c = c+(J) so that 
there is a heteroclinic orbit connecting the saddle point 



(£i,0) to the saddle point (E 3 ,0) with F > 0. Similarly 
c = c_(J) corresponds to a heteroclinic orbit connecting 
(£ 3 ,0) to (£i,0) with F < 0. The functions c±(J) for 
our model are depicted in Fig. 0. Note that they intersect 
when J — J*, given by 

J* = rp 1 „ f 3 v(E) dE, c± = J*, (13) 

-&3 — Hj\ J Ex 



the so-called equal-areas rule pj . 

In the limit 5 <C 1, we can obtain explicit expressions 
for c± (J) and for the corresponding wavefronts as we now 
recall 



The wavefront (moving depletion layer) 
joining (E 3 (J),0) and {E\(J), 0) may be approximated 
by the exact solution of ( J12| ) 



HO = F(0 = -l, 



-{J) = J, 



(14) 



which holds for any value of 5 plus two corner layers of 
width 0(VS) (on the £ scale) at £ = ±(E 3 - E^/2. See 
Ref. |L8) for an explicit calculation. The width of this 
wavefront on the £ length scale is (E 3 — E\) + 0(\/S), 
which yields y — Yi(s) = 0(e) on the large length scale. 
The velocity of this wavefront is c_(J) = J. 

The other wavefronts (moving accumulation layers) 
can be constructed by matched asymptotic expansions 
in the limit S <C 1 and their velocities c + (J) and shapes 
depend on whether J is larger or smaller than J* . These 
wavefronts are composed of a shock wave joining two field 
values £L and E + (at least one of them should be equal 
to Ei, i = 1, 3) plus a tail region which moves rigidly with 
the same velocity as the shock |l(J . The inner structure 
of the shock wave (for very small but not zero S) can be a 
quite complicated triple-deck set of boundary layers . 
Let V(E+, E-) be the velocity of the shock wave given 
by the equal-areas rule [|l0|Jl7|,[lq1 



V(E + ,E_ 



We now have [hOl 



1 



E+ -E. 



v(E) dE. 



(15) 



1. If J G (v m , J*), E + = E 3 (J) whereas E_ is calcu- 
lated as a function of J by imposing the condition 
that the tail region to the left of the shock wave 
moves rigidly with it 



V(E 3 ,E_) 



(16) 



Solving simultaneously ( |15| ) and ([L6|) (with E + = 
E 3 ), we find both E^ and c + = v(E-) > J as 
functions of J. To the left of the shock wave (in 
the tail region) the field satisfies the (approximate) 
boundary value problem 
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[ v (E)-c + (J)]^ = J-v(E) 



E(-oo) = E X {J), 



£<0, 
E(Q)=E_(J). (17) 



2. If J e (J*,vm), E- = E\{J) whereas E + is calcu- 
lated as a function of J by imposing the condition 
that the tail region to the right of the shock wave 
moves rigidly with it 



V(E + ,E 1 )=v{E+). 



(18) 



Solving simultaneously ( |15| ) and (18) (with E- = 
Ei), we find both E + and c + = i>(-E+) < J as 
functions of J. To the right of the shock wave (in 
the tail region) the field satisfies the (approximate) 
boundary value problem 

dE 



[v(E)-c + (J)] — = J-v(E) 



£>0, 

E(0)=E+(J), E(oo)=E 3 (J). (19) 



At J = J* we have £L = E u E + = E 3 , and c + = J*. It 
is not hard to prove that c+(J) is a decreasing function. 
The functions c±(J) are depicted in Fig. |7[ 

Notice that the inner structure of the shock wave has 
width 0(y/S) on the £ scale while the total width of the 
wavefront is 0(1) on the £ scale and 0(e) on the y scale. 
In the limit 5 -C 1, the structure of the wavefronts is 
thus one-sided: the wavefront is a discontinuity preceded 
or followed by a tail region. 

In conclusion, the propagation of a single wavefront can 
be described by giving the position, Y(s), and velocity, 
c±(J(s)), of the front and the values of the electric field 
on its left and right hand sides, Ei(J(s)). When more 
than a single wavefront are co-moving inside the sample, 
the same description applies to each of them. All of the 
magnitudes involved in this description depend on time s 
only through the instantaneous value of the current J(s). 
This fact, together with the fact that the voltage must 
remain constant all the time, can be used to derive a 
simple closed equation describing the time evolution of 
the current density during the propagation stages. As an 
example of this result, let us consider the case a single 
propagating dipole (two fronts) (see Fig. 0b. 1). In this 
case, neglecting transition and boundary layers, we have: 
E(y, s) = Ex(J(s)) for < y < Yi(s) and Y 2 (s) < y < 1 
and E{y,s) = E 3 (J(s)) for Yy > y > Y 2 {s). For the 
voltage we have: 

4> = Bi(J(«)) + [E 3 (J(s)) - 

E 1 (J(s))](Y 2 (s)-Y 1 (s))+0(e). (20) 

By using the fact that the voltage is fixed, we can obtain 
an equation for J(s) by differentiating the bias condition 
( pfj| ) with respect to s. By noting that 



1 



dY x dY 2 



v'(E(J)) ' 
c_(J), 



(21) 



the following simple closed equation for the current is 
obtained 



^=A(J) [c+(J)-c-(J)] 



A(J) = 



(E 3 — Ei) 2 



>0, 



(22) 
(23) 



where v[ = v'(Ei(J)) (i — 1,3). This equation for J 
holds as long as the electric field profile consist of a single 
propagating dipole. A simple analysis of Eq.(|2^) demon- 
strates that J tends to J* (for which c+(J) = c_(J)) 
exponentially fast, starting from a certain value of the 
current J(0). This result explains, in particular, why the 
solitary wave responsible for the Gunn effect moves at 
approximately constant velocity c± = J* and constant 
current far from the contacts. 

A similar procedure can be applied to derive the corre- 
sponding closed equation for the current when different 
types of propagating domains are present. In general we 
will have n + fronts moving at velocity c + and n_ moving 
at c_, with |n+ — n_| = 0, 1. For these cases the current 
evolves following the equation 



dJ 
ds 



= A(J) [n+c+(J) - n_c_(J)] . 



(24) 



As before, this equation holds as long as the propagating 
profile consists of n + fronts propagating at c + and n_ 
at c_. Two typical cases can be considered in analyzing 
Eq. (24): (i) cither n + or n_ are equal to zero and (ii) 
n± 7^ 0, both numbers are different from zero. In the 
first case, we have either 



or 



f s =AiJ)c +{ J), 



d J. = -A(J)c A J), 



n+ 



1 , n_ 



Q,n_ = 1 



(25) 



(26) 



and the current will either increase, Eq. (25), or decrease, 
Eq. (|26|), with time. In the second situation, the cur- 
rent follows the general relation Eq. (^) with the cor- 
responding values of n±: the current evolves towards 
the fixed point J — J n+ ,„_, satisfying n+c+(J n+ . n _) — 
rt_c_( J n+jn _ ) — 0, when such a point exists. A partic- 
ular example of this behavior has been explained above 
for Ji,i = J*. 
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These results describe the quasistationary propagation is given by 
of fronts and the corresponding time evolution of the cur- 



rent density during these stages. We will use them to ( — + P J [J^ — h'(a) — ah] = —7/1, (29) 

interpret the results of the numerical simulations. \do~ J 



B. Generation of fronts 



As mentioned before, for our b.c. an appreciable part 
of a period of the current oscillation may be spent gen- 
crating new fronts non-adiabatically. Fronts of dipolc 
domains are generated at the cathode whereas monopole 
wavefronts appear somewhere in the middle of the sam- 
ple. In what follows, we will focus on the description of 
dipole domains. Monopole wavefronts have been recently 
described in detail elsewhere JTT| . 

A simple rule concerning generation of dipole domains 
can be formulated: when the current density, J(s), 
crosses the maximum [resp. minimum] of the contact 
curve, j C M [resp. jcm], a front moving with speed c_ [resp. 
c+], starts being formed when its generation is compatible 
with the field value at the bulk after the injecting contact. 

Let us now show why the previous rule holds. Sup- 
pose that the current reaches adiabatically in the slow 
time scale s one of the critical values mentioned above, 
let us say j C M- Then the field at x = 0, given by Eq. 
(^), can no longer be quasistationary, the injecting layer 
becomes unstable and it starts shedding a new solitary 
wave. Let now s\ be the earliest time at which J = j C M, 
and the boundary field Eq = Em- After this time, the 
disturbances J — j C M = 0(e) cause Eq to evolve to the 
third branch of j c (E) on a time of order s — Si = O(e). 
The field in the injecting layer, in turn, increases until a 
moving wavefront moving with speed c_ is formed. To 
describe this process, we adapt ideas developed for the 
analysis of the trap-dominated Gunn effect in p-Ge [ 12] ] 
to the present situation. An important difference is that 
the sharp increase in the contact field Eq helps creating 
the new wavefront but a new pulse is not shed in the fast 
time scale a = (s — s±)/e: after the wavefront is created, 
the contact field varies on the third branch of j c (E) and 
the current has to decrease slowly [according to (|3), on 
the s time scale] until Eq can jump back to values on the 
first branch of j c (E) when J = j cm . 

To leading order, the field in the injecting layer solves 
the following semi-infinite problem whose derivation can 
be found in Appendix H: 



dE (BE 
— +v(E) — 
oa \ ox 



1 



J, 



x > 0, —00 < a < +00 , 

E(0,a) = Eq(o-), 



(27) 



(28) 



where h(a) is 

poo 

h(a) = (E 3 -E 1 )c+a- [E{x, a) - E x ] dx 

Jo 

-f [£(£)- SiO'cmM 

J —OO 

) 

[E(0-E 3 (j cM )}d^, (30) 



and formulas for the pos itive constants a, (3 and 7 may 
be found below, in Eqs. ( B1C ). The function h(a) is the 
area lost due to the motion of the old front during the 
time a minus the instantaneous excess area under the in- 
jecting layer minus the constant excess area under the old 
wavefront at Y = Yx(sx) [that is, under the heteroclinic 
orbit connecting (Ei(j cM ),Q) and (E 3 (j cM ),0)}. 

As a — > — 00, we have to impose the following matching 
condition on an appropriate overlap domain: 



E(x,a) -E stat (x; J(s)) < 1, 
as <y — > —00, s — > si — . 



(31) 



where Eq solves (||) with J = j c m + eJ^ 1 ^(cr), and J^\a) 



Here E sta t{x; J(s)) is the quasistationary injecting layer 
solution of ( |l2| ) with c = such that E sta t(0', J{s)) sat- 
isfies (§) and E stat (oo; J(s)) = Ei(J(s)) for s < si, 
>^( s i) — icM- Notice that the term eJ^'(cr) is needed 
so as to avoid that the solution of this problem stay in- 
definitely in the quasistationary field E stat (x; J(s)). 

The solution of the previous semi-infinite problem re- 
veals the growth of the field at the contact and in- 
side the injecting layer until: (i) Eq becomes E c3 (j c m) 
[E cl (J) < E c2 (J) < E c3 (J) are the three possible solu- 
tions of J = j c (E)], (ii) E(x,a) increases to E 3 (j c m) as x 
increases from x — and then (iii) it has the structure of 
a wavefront connecting E 3 (j C M) to £a(j C A/)- This wave- 
front advances with velocity c_(j c m)- The formation of 
a front when J crosses j cm can be explained similarly. 



C. Putting the pieces together 



In the two previous subsections, we have presented the 
basic features of the asymptotic description of the Gunn 
effect, namely, quasistationary front propagation and the 
generation of new fronts. Now, we are in a position to 
put all the pieces together, and describe a full period 
of current oscillation. We shall distinguish three cases: 
high-field dipoles, low-field dipoles and monopoles. 



G 



1. Dynamics of high-filed dipoles 



High-field dipolc domains have been observed to ap- 
pear when v m < j cm and J* < jcU < % ■ Then depend- 
ing on the value of j cm with respect to J* , and of j C M 
with respect to ,P = J2.i1 different situations may occur. 

Let start considering the case v m < j cm < J* < j C M < 
Jt < vm- This situation corresponds to the propaga- 
tion of high-field dipole domains as shown in Fig. H. Wc 
will assume the initial electric field configuration to cor- 
respond to a single propagating high field domain (Fig. 
|l|b.l). This configuration corresponds to n + — 1 and 
ri_ = 1, and hence the current satisfies Eq. (p2|), evolv- 
ing from an initial value J(0) € (jcm, J*), towards the 
fixed point J* = Ji,i- After a certain time, the wavefront 
located at Y2 reaches the end of the sample and disap- 
pears, producing an abrupt change in the time deriva- 
tive of the current. We have a new stage with n + = 1 
and n_ = 0, Fig. 0b. 2, governed by Eq. (|25|). Its solu- 
tion increases until it surpasses the value j c M- At this 
point the injecting layer becomes unstable and it starts 
shedding a new front. The formation dynamics has been 
explained in detail in the previous section. Then a new 
slowly varying stage begins. There are two leftover wave- 
fronts, the old one located at y = Yi(s) [which advances 
toward y = 1 with speed c+(J)], and a new one located 
at y = Yi(s), moving with speed c_(J), and leaving be- 
hind a quasistationary field Es(J), see Fig. |l|b.3. Again, 
the field configuration corresponds to n + = n_ = 1, with 
J(s) following (p2]), and decreasing exponentially fast to- 
wards J*. Before J reaches J*, the front located at Y\ 
reaches the end of the sample and disappears, thereby 
producing a new abrupt change in the time derivative of 
the current density. Then only the recently formed front 
[located at ^4(5)] is present on the sample (Fig. [l]b. 4), 
which corresponds to n + — and n_ = 1. J(s) de- 
creases according to Eq. (^6|) until the minimum of the 
contact curve, J = jcm, is reached. After that, a front 
moving with speed c+ is generated. The charge dipole 
wave thus created evolves adiabatically, and the current 
density is described again by Eq. (§§), Fig.03.1. We have 
come back to the initial situation, and one period of the 
current oscillation has been completed. 

It is worth noting that by means of this analysis some 
of the most relevant features of the current oscillations, 
as the maximum and minimum currents, Jmax and Jrnim 
have been identified with quantities related to the con- 
tact parameters, J max ~ jcM and J mi „ ps j cm , opening 
then the possibility of determining the values of contact 
parameters from the analysis of the current oscillations. 

Other situations involving dipole domains have been 
numerically identified. They are described by the same 
type of analysis. In what follows only the more relevant 
features of these situations will be considered. Let us as- 



sume now v m < J* < j cm < j C M < J' < vm- This case 
corresponds again to high-field solitary waves (dipole do- 
mains, see Fig. ||), but the current oscillations have a dif- 
ferent shape. The three first stages of the oscillation, Fig. 

.1,2,3, correspond to the propagation of a single dipole, 
annihilation of a front and generation of new one. They 
are similar to the stages described above. Now, however 
after the new front has been formed, Fig. |]b.3, and the 
current is decreasing towards J*, J reaches the critical 
current, j cm (because J* < j cm )- A new front moving 
with speed c+ is then created. After the formation pro- 
cess has finished, we have a configuration with two fronts 
moving at c+ and one at c_, that is, with n+ = 2 and 
n_ = 1 (Fig. |]b.4). Hence the current satisfies [see Eq. 

©], 

^=A(J)[2c + (J)-c_(J)], (32) 

and it starts increasing again, trying to reach J\ t % = JL 
Before this value may be attained, the old front located 
at Y\, arrives at the receiving contact and disappears. 
We obtain again Eqs. (53)- (pa) and recover the initial 
situation. Thus a full period of the Gunn oscillation is 
again completed. 

There are other possible situations for propagating 
high field domains, but we have not found them in our 
numerical simulations with the curve v(E) considered in 
the present work. For example in the second case we have 
described above, v m < J* < j cm < j c M < < vm, af- 
ter the formation of the new front, Fig. Hb.4, the current 
density could have reached the value J< before the old 
front located at Y\ had arrived at the receiving contact. 
In such situation, a new front moving at c_ could be 
formed, giving rise to an electric field configuration with 
n + = n_ = 2. In this configuration the current would 
decrease again to J2.2 = Ji,i — J*- We would then have 
a more complicated situation with two pulses and the 
dying wavefront simultaneously present in the sample. 
Proceeding in a similar way, we could therefore observe 
simultaneous coexistence of several pulses for appropriate 
ranges of parameters. This situation (multiple shedding 
of high field domains) has been numerically observed in 
other models 



2. Dynamics of low-field dipoles 

Low-field dipole domains appear when v m < j cm < 
jcM < J* (%•)• Depending on the value of the applied 
bias, single (high bias values) or multiple (low bias val- 
ues) propagating low field domains are obtained [see Figs. 

Let us start considering the case of single propagating 
low field dipole domains, Fig. 0. We consider an initial 
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field configuration having a low field domain far from the 
contacts. Then there is a wavefront located at Y\, mov- 
ing with speed c_, and a wavefront located at Y2 > Y\, 
moving with speed c+, that is, n + = n_ = 1, Fig. |3|b.l , 
and J(0) € (j m i J*)- The current increases towards J* 
according to Eq. (|22|), until Y2 — 1. Then J starts satis- 
fying (26) (corresponding to n+ = and nt = 1, see Fig. 
gb.2), and it decreases until the value j cm is reached. 
After this occurs, a new front (moving at speed c + ) is 
created, while the old wavefront Yi < 1. We have a stage 
described by Eq. ( p2| ) with ?i + = n_ = 1, see Fig. ||b.3, 
during which J increases past j C M {jcM < J*)- Then a 
new front moving with speed c_ starts being created. At 
the same time, the old front located at Yi reaches the end 
of the sample, giving rise to a complex stage in which the 
non-stationary effects can not be neglected Fig. Hb.4. Af- 
ter that stage, we recover the first situation Fig. kjb.l, and 
a complete period of the oscillation has been described. 

Let us now consider an example of multiple propagat- 
ing low-field dipole domains (Fig. ^). As we mentioned 
above, they appear for small bias values. The asymptotic 
analysis of this case is more complicated because there 
are new stages (fusion of two wavefronts inside the sam- 
ple) whose detailed description is outside the scope of 
this paper. Let us start with the configuration shown in 
Fig. |b.l, for which two low field domains coexist. This 
configuration corresponds to n + = n_ = 2. Following 
Eq. (E4T) , the current will evolve according to 



dJ 
ds 



A{J) [2c + (J)-2c_(J)], 



(33) 



thereby increasing towards Ji,2 = J* ■ Before this value 
can be reached, the front located at Y4 reaches the end 
of the sample. The resulting configuration (Fig. ^b.2) 
has now n + = 1 and n_ = 2, and the current decreases 
according to 



^=A(J) [c+(J)-2c_(J)]. 



(34) 



If the curve v (E) were such that the fixed point of Eq. 
(ID) existed, the current would tend to such a value. For 
our choice of v(E), this fixed point does not exist (see 
Fig. [7]), and therefore the current decreases all the time. 
During this process j cm will be crossed. According to 
our previous considerations, a new front moving at speed 
c+ should then be formed. However, this does not oc- 
cur because melting of the fronts Y2 and Y3 starts after 
Y4, reaches y = 1. This melting process seems to in- 
hibit the creation of new wavefronts until it is completed, 
which happens for J < j cm . Then a new front moving 
at speed c + is rapidly created and we have a stage, with 
?!+ = n- = 1, Fig. |b.3, described by Eq. (§|). The 
current increases until J = j C M, at which time a front 
moving at speed c_ appears at y = 0, Fig. ||b.4. Then 
the current decreases following (|I|) until J = j cm . An- 
other front moving at speed c+ is then formed. We have 



now 7j + = n_ = 2, Fig. Qb.5 and J increases according 
to (p2|). When J > j C M, a front moving with speed c_ is 
created at y = 0, so that n+ = 2, n_ = 3, see Fig. ||b.6. 
Now the current decreases towards J2,3 according to the 
equation 



d£ 
ds 



A{J) [2 C+ (J)-3c_(J)]. 



(35) 



This stage lasts until Y\ = 1, at which time we are back at 
the first stage, Fig. ||b.l, having completed a full period 
of the oscillation. More complicated examples of multi- 
ple low field domains exist and they could be described 
similarly. 



3. Dynamics of charge monopoles 



When j cM > v M , or when j cm < v m &ndj sat > v m , the 
Gunn effect is mediated by wavefronts which are charge 
monopoles. The case j C M > vm, corresponds to mov- 
ing accumulation layers, Fig. |^, while moving depletion 
layers are observed for j cm < v m and j sat > v m , Fig||. 
H. Kroemer fl5|| discovered numerically the Gunn effect 
mediated by charge accumulation monopoles for bound- 
ary conditions corresponding to the control characteris- 
tics, while its asymptotic analysis (for B« 1) was per- 
formed many years later [Tc| ]. Recently a Gunn effect 
mediated by charge accumulation monopoles has been 
used to describe self-sustained oscillations of the current 
in GaAs/AlAs superlattices in the limit of weakly 
coupled, low doped, long superlattices |LjJ . 

The complete asymptotic analysis found in Ref. |Tl| ] 
can be used to describe the present situation after a few 
simple changes are made. First of all the equal-areas 
rule in the case of the superlattice refers to l/v(E), not 
v(E) p2fl . The second important change is that of the 
injecting boundary condition in the description of the 
birth of a new monopole: instead of a rigid Neumann 
boundary condition dE(0,t)/dx = c we should use (Q). 
These two changes can be implemented without difficulty 
and the details will be omitted. An important difference 
with the Gunn effect mediated by dipole solitary waves 
is that the amplitude of the current oscillations is larger 
(its largest value is approximately vm ~ v m for the charge 
accumulation monopoles and jcM — jcm for the dipoles). 
The monopoles "probe" the full region of negative slope 
of v{E) while the dipoles "probe" a smaller region. 



IV. CONCLUSIONS 



We have performed an asymptotic and numerical anal- 
ysis of the Gunn effect in n-GaAs under general boundary 
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conditions for metal-semiconductor contacts. We have 
shown that the Gunn effect is mediated by (i) mov- 
ing depletion charge monopoles, (ii) moving accumu- 
lation charge monopoles, (hi) high-field dipole solitary 
waves or (iv) low-field dipole solitary waves, according 
to whether the critical contact currents j c m = 3c{Em), 
jcm = jc(E rn ), j^ at are (i) j cm < v m and j^ at > v m , (ii) 
jcM > v M , (ih) v m < jcm and J* < j cM < v M , or (iv) 
v m < jcm < jcM < J*, in dimensionless units [v m and 
Vm are the minimum and maximum values of the elec- 
tron drift velocity v(E), E > 0, and J* is the solution of 
c + (J) = c_(J)]. Some of these results are well-known for 
boundary conditions given by Kroemer's control charac- 
teristic. In addition we have shown that there are new 
instability mechanisms consisting of multiple generation 
of (low-field) charge dipole solitary waves in the region 
near the injecting contact if j C M is close enough to J*, 
and the dimensionless length is large enough. In each 
case we have been able to describe in detail the shape of 
the current oscillation, identifying some of its main fea- 
tures (maxima, minima, plateaus, . . . ) with critical cur- 
rents appearing in our model, among them the contact 
currents j cm and j C M- This result opens the possibility 
of determining contact parameters by means of the anal- 
ysis of the shape of the current oscillation. Our results 
might be of use in the analysis of self-sustained current 
oscillations in weakly coupled n-doped superlattices [ElJ , 
once the role of contacts and the boundary conditions 
they generate are understood in these systems. Beside 
this, our results are not restricted to the particular model 
of the Gunn effect studied here, but seem to hold for a 
general class of models, supporting the idea of the "uni- 
versality" of the Gunn effect |i~3|] . 



This can be seen from the method of characteristics ap- 
plied to (|l|) with 5 = and the boundary condition (|) 
The characteristic equations are 



§ -'-<*>• 

!-•<*>■ 



to be solved with the conditions 

E{t-t) =E (t), 
x(t; t) = 0. 



(A2) 
(A3) 



(A4) 
(A5) 



Clearly, 



dEo 
dr 



dE. . dE. . 
— (t;t) + —(t;t) 



= J -v(E q ) + —(t;t). 



(A6) 



The last term in this equation can be obtained from the 
boundary condition (||) as follows. From the solution of 
(JO) and (|A|), we obtain ff(r;r) = -v(E (t)). Then 
the boundary condition (Eh yields: 



dE , s dE , Bx, x 
-v(E q (t)) a (i - J + — ) . 



(A7) 



Insertion of (A7) into (|Aq) yields (Al) (with r = t) 
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To obtain the equations governing the shedding stages, 
we need to keep terms of order e in our calculations, as 
indicated in Ref. [fT2[. The outer (bulk) expansion is 



E = E?\a)+eEY > {o)+0{e% 
J = J ( -°\a)+eJ^(a)+0(e 2 ), 



(Bl) 
(B2) 



APPENDIX A: NONSTATIONARY FIELD AT 
THE INJECTING CONTACT 



where i = 1 if < y < Yi(si) and i = 3 if Yi(si) < y < 1. 
Then K (0) and solve 



We shall now prove that the field Eq (t) at the contact 
x = obeys the equation: 



dE, 



(0) 



da 



+ v(E^) = j(°\a), 



E 



(0) 



{E^-Emi-Y^)}. 



,(0) 



(B3) 
(B4) 



- W = J-] C (E ), 



(Al) 

We can obtain easily E$ and from these equations: 







E. 



(0) 



J<°>(<t) = v(E[ 0) )Y 1 ( Sl ) + v(E^)[l - 



(0)x 



(B5) 
(B6) 



Now £^°' ) > follows from (B3) and these two last equa- 
tions. These equations should be solved with the match- 
ing conditions that — ► Ei(j C M) {i = 1,3) and 

J(0) 

— ► j C M as a — > — oo. These problems have the solu- 
tions J = and El = Ei(j c M)- Thus we expect that 
the current density does not depart substantially from 
j cM as e — > 0. 

The O(e) corrections E"} 1 ^ (cr) obey the equations: 
Q t E^ ee ^ +v ' {El )E[ 1] = J«(a), (B7) 



which immediately follow from (JsJ) . To find we pro- 
ceed as follows. First of all, we write the bias condition 
including terms of order e in the approximation of 
the electric field and the current density: 

<t> = E x {j cM ) + [E 3 {j cM ) - Ei{j eM )] [1 - Y 1 (s 1 )} 
+e {e[ 1} + - B«] [1 - Yi(si)] - [S 3 (i cM ) 

/•oo 

--E'i(jcM)]c+cr + / cr) - Ei(j cM )] dx 



( [E(0-E 1 (j cM )}d4 

J —oo 



+ J [E(£)-E 3 (j cM )]di}+0(e*). (B8) 



Here E(£) is the field inside the wavefront at Y\ and 
E(x, a) is the field in the injecting layer. Since <f> — 
Ei(J c m) + [E 3 (j cM ) - E^jcu)} [1 - Yi(«i)], we obtain 



All functions of J in these equations are calculated at 
J = jcM- Solving this equation we obtain the following 
function 



J (1) (fj) = h'(a) + ah(a) - 7 



h{t)dt. (Bll) 



The only missing function is the field at the inject- 
ing boundary layer. This field profile is the solution of 
the semi-infinite problem (£7 ) to (|3l|), with J (a; e) = 
jcAf + eJ^^a) given by (Bll). We can write ( Bll ) in 
another form that suggests a more transparent interpre- 
tation: 

jW{a) = ~{J'( Sl )(a-a )-l'(a) 

- J™ [aI'{o - a') + v[v' 3 I{a - a')} e^"' da'\ , (B12) 
f°° 

1(a) = - / [E(x, cr) - EA dx, (B13) 
a Jo 

+{[E3(jcM)~E 1 ( JcM )}c + }- 1 
xf[ [E(g) - Ex{j cM )]d£, 

K.J — OO 

+ J^°[E(0-E 3 (j cM )}A. (B14) 



CT = 



The terms on the right side of ( B12 ) clearly display the 
balance between the area lost by the motion of the old 
wavefront at Y\ (s) and the excess area under the injecting 
layer. J (1 V) 

increases linearly with a until the excess 
area under the injecting layer increases with a at least 
linearly. 



E[V + [l-Y 1 (s 1 )]=h(a), 



(B9) 



where h(a) is given by (J3 1 

We now obtain the first order differential equation J29|) 
for jW(cr) by applying the operator 6163 to both sides 
of (p9) and then using (B7): 



d_ 



+ 0) [J (1) - h'(a) - ah] = -7/1, 



where the coefficients a, and 7 are 

a = v[Y 1 (s 1 ) + [l-Y 1 (s 1 )]v 3 , 
= v' 3 Y 1 (s 1 ) + [l-Y 1 (s 1 )]v' 1 , 

7 = ri( S i)[i-ri( S i)]K-« 3 ) 2 . 
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FIG. 1. Gunn effect mediated by dipole solitary waves, 
(a) Dimensionless current density J(s). Parameter values 
are L — 800, io = 0.3, o?o = 3.6 and 4> — 1.5, for which 
v m < jcm < J* < jcM < vm- The minimum and maximum 
values of J(s) corresponds to j cm = 0.23 and j C M = 0.31, re- 
spectively whereas the plateau at intermediate values of J(s) 
corresponds to the solution of c+(J) = c_(J), J* = 0.28. (b) 
The corresponding electric field profiles E(y, s) evaluated at 
the times marked in Part (a) of this figure. 
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FIG. 2. Gunn effect mediated by dipole solitary waves (a) 
Dimensionless current density J(s). Parameter values are 
L = 1000, i = 0.45, a = 30 and = 2.4, for which 
» m < J* < j cm < j C M < vm- The maximum value of J(s) 
corresponds to j C M = 0.46, whereas the plateau at bottom of 
J(s) corresponds to the solution of c+(J) = c_(J), J* = 0.28. 
The value of j cm is 0.41. (b) The corresponding electric field 
profiles E(y, s) evaluated at the times marked in Part (a) of 
this figure. 
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FIG. 3. Gunn effect mediated by low-field dipole solitary 
waves, (a) Dimensionless current density J(s). Parameter 
values are L — 1500, io = 0.24, ao = 30 and <f> = 2.3, for 
which v m < j cm < jcM < J* < vm- where J* = 0.28 cor- 
responds to the solution of c + (J) = c_(J). The values of 
jcM = 0.26 and j C m = 0.23. (b) The corresponding electric 
field profiles E(y, s) evaluated at the times marked in Part 
(a) of this figure. 
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FIG. 4. Multiple shedding of low-field domains. The pa- 
rameter values are L — 800, io = 0.24, a = 30 and 
cj> — 1.5. The critical currents are J* = 0.28, j C M = 0.26 
and jcm = 0.23. (a) Dimensionless current density J(s). (b) 
The corresponding electric field profiles E(y, s) evaluated at 
the times marked in Part (a) of this figure. Two pulses are 
formed during each period. The secondly shed pulse reaches 
and overtakes the first one. 
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FIG. 6. Gunn effect mediated by monopole wavefronts 
(moving depletion charge monopoles) . (a) Dimensionless cur- 
rent density J(s). Parameter values are L — 800, io = 0.135, 
ao — 0.8, and cj> — 1.5, for which j C M > vm- (b) The corre- 
sponding electric field profiles E(y, s) evaluated at the times 
marked in Part (a) of this figure. 
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FIG. 5. Gunn effect mediated by monopole wavefronts 
(moving accumulation charge monopoles). (a) Dimension- 
less current density J(s). Parameter values are L = 800, 
io = 1.35, ao = 0.8, and (j> — 2, for which j cm < v m and 
jo at > v m , where jQ at corresponds to the saturation current 
jo at = a -1 + io- The maximum and minimum values of the 
oscillation correspond to the maximum and minimum values 
of the v(E) curve, vm = 0.58 and v m — 0.18, respectively, 
(b) The corresponding electric field profiles E(y, s) evaluated 
at the times marked in Part (a) of this figure. 
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FIG. 7. Velocities c+ and c_ of the heteroclinic wavefronts 
as functions of the current density. Notice that the lines 
intersect at J* = 0.28. We have also marked the currents 
— J2,i = 0.53 at which 2c + = c_ and J-2,3 = 0.20 at which 
2c+ = 3c_. 
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